function [model] = update_alpha_02(model)

%global model;

N              = model.N;
term           = (model.alpha_2 - model.lambda*model.alpha_1*model.eta');
model.alpha_02 = [sum(term,1)/N];

end
